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Abstract. 

We study the finite-size effects for the thermal QCD Deconfinement Phase Transition (DPT), and 
use a numerical finite size scaling analysis to extract the scaling exponents characterizing its scaling 
behavior when approaching the thermodynamic limit (V — > oo). For this, we use a simple model 
of coexistence of hadronic gas and color-singlet Quark Gluon Plasma (QGP) phases in a finite 
volume. The Color-Singlet Partition Function (CSPF) of the QGP cannot be exactly calculated 
and is usually derived within the saddle point approximation. When we try to do calculations with 
such an approximate CSPF, a problem arises in the limit of small temperatures and/or volumes 
' (VT 3 << l), requiring then additional approximations if we want to carry out calculations. We 

, propose in this work a new method for an accurate calculation of any quantity of the finite system, 

^"*^ ■ without explicitly calculating the CSPF itself and without any approximation. By probing the 

behavior of some useful thermodynamic response functions on the hole range of temperature, it 
ON , turns out that in a finite size system, all singularities in the thermodynamic limit are smeared out 

and the transition point is shifted away. A numerical finite size scaling analysis of the obtained data 
allows us to determine the scaling exponents of the QCD DPT. Our results expressing the equality 
between their values and the space dimensionality is a consequence of the singularity characterizing 
a first order phase transition and agree very well with the predictions of other FSS theoretical 
approaches and with the results of both lattice QCD and Monte Carlo models calculations. 



INTRODUCTION 



o 

It is generally believed that at sufficiently high temperatures and/or densities a new phase of matter called 
j^H. the Quark Gluon Plasma (QGP) can be created. This is logically a consequence of the quark-parton level of the 
matter structure and of the dynamics of strong interactions described by the Quantum Chromodynamics (QCD) 
theory. Its existence, however, has been partly supported by lattice QCD calculations and the Cosmological 
Standard Model predictions. The only available way it seems to study experimentally the QCD Deconfinement 
Phase Transition (DPT) is to try to create in Ultra-Relativistic Heavy Ion Collisions (URHIC) conditions similar 
to those in the early moments of the Universe, right after the Big Bang. If ever the QGP is created in URHIC, 
the volume within which the eventual formation would take place would certainly be finite. Also, in lattice 
QCD studies, the scale of the lattice space volume is finite. This motivates the study of the Finite Size Effects 
(FSE) on the expected DPT from a Hadronic Gas (HG) phase to a QGP phase. These effects are certainly 
important since statistical fluctuations in a finite volume may hinder a sharp transition between the two phases. 
Phase transitions are known to be infinitely sharp, signaled by some singularities, only in the thermodynamic 
limit yj. In general, FSE lead to a rounding of these singularities as pointed out in |2j,|3|. However, even in a 
such situation, it is possible to obtain informations on the critical behavior. Large but finite systems show a 
universal behavior called "Finite-Size Scaling" (FSS), allowing the put of all the physical systems undergoing 
a phase transition in a certain number of universality classes. The systems in a given universality class display 
the same critical behavior, meaning that certain dimensionless quantities have the same values for all these 
systems. Critical exponents are an example of these universal quantities. 

In the present work, we study the finite size behavior for the thermally driven QCD DPT. For this purpose, 
we consider a simple Phase Coexistence Model (PCM) used in [4(, in which the mixed phase system has a 
finite volume: V = Vhg + Vqgp- The parameter t) representing the fraction of volume occupied by the HG: 
Vhg = fyV, is then defined and can be considered as an order parameter for the QCD DPT. Assuming non- 
interacting phases (separability of the energy spectra of the two phases), and using the PCM, the total partition 
function of the system can be written as: 

Z(fj) = Z QGP ty)Z HG ty)Z Vac $), (1) 
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where: 

Z Vac (T, f)) = exp(-BV QGP /T) = exp(-W(l - f))/T), (2) 

accounts for the confinement of the color charge by the real vacuum pressure exerted on the perturbative vacuum 
(B) of the bag models. 

The mean value of any thermodynamic quantity of the system (A(T, [i, V)}, as defined in can then be 
calculated by: 

fA{t),T,n,v)z(t))di) 

(A(T,n,V)) = * j , (3) 

o 

where A(t), T, /x, V) is the total thermodynamic quantity in the state f), given in the case of an extensive quantity 
by: 

AQ),T,»,V) = Ahg{T,ii,W) + Aqgp{T,h, (1 - f))V), (4) 
and in the case of an intensive quantity, by: 

A(t), T, (x, V) = i)A HG (T, fi, i)V) + (1 - t))A QGP (T, fi, (1 - t))V), (5) 

with Aqgp and Aug the contributions relative to the individual QGP and HG phases, respectively. 

We think that the PCM used in our present work is not a characteristic signature in a finite volume which 
anticipates the behavior of a first order phase transition in the TL. It is well known that finite size effects lead 
to a non obvious order of the phase transition. A characteristic feature of a truly first order transition in the 
TL can then appear in the second order transition when the volume of the system is finite and vice versa, and 
this is considered as a finite size artifact. For example, the phenomenon of coexistence of phases can appear in 
a truly second order phase transition when the volume of the system is finite, and we can also have a divergent 
correlation length in a truly first order phase transition Certainly, these non-characteristic features must 
disappear when approaching the TL. 

For the HG phase, we have considered just pionic degrees of freedom, and the partition function for such a 
system is simply given by: 

Z HG =e^ T3yHG . (6) 

For the QGP phase, we have considered a free gas of quarks and gluons with the exact color-singlctness 
requirement. The color-singlet partition function of the QGP derived using the group theoretical projection 
technique formulated by Turko and Redlich |g] can not be exactly calculated and is usually calculated within 
the saddle point approximation in the limit VqgpT 3 >> 1 as in 0, 0, El ■ ^ turns out that the use of the 
obtained approximated partition function for the calculation of a mean value within the definition (3), has 
as a consequence the absence of the DPT ^(J- This is due to the fact that the approximation used for the 
calculation of the color-singlet partition function breaks down at VqgpT 3 << 1, and this limit is attained in 
our case. This has been emphasized in further works in which additional approximations have been used to 
carry out calculations, as in 0,0- We propose in the following a new method which allows us to accurately 
calculate physical quantities describing well the QCD DPT at finite volumes, avoiding then the problem arising 
at VqgpT 3 << 1 and without any approximation. We proceed by using the exact definition of the color-singlet 
partition function for Zqgp, without explicitly calculating it, in the definition (3) of the mean value of a physical 
quantity. A first analytical step in the calculation of the mean value is then achieved, and an expression with 
integral coefficients is obtained. The double integrals are then carried out with a suitable numerical method at 
each value of temperature and volume, and the behavior of the physical quantity of the finite system can so be 
obtained on the hole range of temperature, for various volumes, without any restriction. Afterwards, scaling 
critical exponents characterizing the scaling behavior of some quantities are determined using a numerical FSS 
analysis. 
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2. EXACT COLOR-SINGLET PARTITION FUNCTION OF THE QGP 



The exact partition function for a color-singlet QGP contained in a volume Vqgp, at temperature T and 
quark chemical potential /i, is determined by 



+ 7T +7T 

■^qgp(T, Vqgp, m) = ^ / / d (f ) d (f ) M ^ ^( T > W V, tf), (7) 

— 7T — 7T 

M((p : ifc) is the weight function (Haar measure) given by: 

M(^)= (sin(i(V; + f))sin(f) S in(i(^-f))) 2 , (8) 
and Z the generating function defined by: 

Z(T, V QG p, w ^) = Tr [exp (-/? (ff - // (X - N^j ) + ipl 3 + tyy 8 )] , (9) 

where (3 = — (with the units chosen as: kg = ft = c = 1), iJo is the free quark-gluon Hamiltonian, N q (j^ijj 

denotes the (anti-) quark number operator, and I3 and Y% are the color "isospin" and "hypercharge" operators 
respectively. Its final expression, in the massless limit, can be put on the form: 

Zqgp(T,Vq G p^) = ^ J * j "rf^^^e^^^.*), (10) 



with: 



^^'T )= 12 ( 30 dQ+ 15 dG)+ 12T D - 1+ ^ 1 

q=r,b,g ^ \ / 



7T 2 <^G / (Q!g - 7r) 2 



^g(^-!) 2 . (ID 

c£q = 2iVy and do = 2 being the degeneracy factors of quarks and gluons respectively, a q (q = r, b, g) the angles 
determined by the eigenvalues of the color charge operators in eq. (9): 

OL r = j + 3 , OL g = — £ + Ckb = 3 , (12) 

and a g (g = 1, ...,4) being: 

ai = a r — a g , a.2 = ct g — a^, = ab — a r , = 0. (13) 



3. FINITE-SIZE EFFECTS 



To study the effects of volume finiteness on the thermal QCD DPT within the PCM, we'll examine in 
the following the behavior of some thermodynamic quantities of the system with temperature, at a vanishing 
chemical potential (/i = 0), considering the two lightest quarks u and d [Nf — 2), and using the common value 
B 1 / 4 = UbMeV for the bag constant. 

The first quantity of interest for our study is the order parameter, which can be simply in this case considered 
as represented by the hadronic volume fraction. According to (3), its mean value is expressed as: 

+ J + f d<pdiPM(ip, V) J qe^^ T - y >dq 

< h(T, V) >= 1 - °- , (14) 

/ f dtpdi/)M(<p, ip) J e^(v^;T,v) d([ 

— IT —7T 
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with: (</?, V 1 ! T, V) — (g fJ _=o( l P,'<P) ~ fo ~ T*) ^^ 3 - After integration on the q variable, the order parameter 



can be written as: 



<mv) >= £ °i+^»- £ " > (is) 

J^m — J^ii 

where the general form of the integral terms appearing in this expression is: 

(16) 

— 7T — 7T 

These integrals are then carried out using a suitable numerical method at each fixed temperature and volume. 
The second quantity of interest is the energy density e(T, V), whose mean value is related to < f)(T, V) > by: 

< e{T, V) >= e HG + (B — e HG ) (1- < J,(T, V) >) - ZT i Ll1 ~ Ll2 + L ° 2 , (17) 

-^01 — Ml 



where the new integrals on ip and if) , noted L nm , are given by: 



+^ r e 9\( v ,i,:T,v)\ n 



— TT J —TT 



and: e HG = j^T 4 . 

We can also examine the effects of volume finiteness by illustrating two more quantities representing the first 
derivatives of the two previous ones, i. e., the susceptibility x defined as: 

, (T , F) = TO>, (19) 

and the specific heat density c (T, V) defined as: 

c(T,V) = 2™ (20) 

Our results for the variations of these four response functions with temperature at different system sizes are 
presented in the plots on Figs. (1) and (2) and show a pronounced size dependence over almost the entire 

< £ > 

temperature range. The curves show that in the limit of an infinite volume, both (f)) and — — exhibit a 
finite sharp discontinuity, which is related to the latent heat of the DPT, at a bulk transition temperature 

T c (oo) ~ 104.35Afey, reflecting the first order character of the transition. The quantity is traditionally 

interpreted as a measure of the number of effective degrees of freedom; the temperature increase causes then 
a "melting" of the constituent degrees of freedom "frozen" in the hadronic state, making the energy density 
attain its plasma value. This finite discontinuity can be mathematically described by a step function, which 
transforms to a <5-function in \ an d c. When the volume decreases, the four quantities vary continuously such 
that the finite sharp jump is rounded off and the <5-peaks are smeared out into finite peaks over a range of 
temperature 8T (V). This is due to the finite probability of presence of the QGP phase below T c and of the 
hadron phase above T c , induced by the considerable thermodynamical fluctuations. 

Another feature which can be noted is that the maxima of these rounded peaks occur at effective transition 
temperatures T c (V) shifted away from the bulk transition temperature for infinite volume T c (oo). This shift 
is a consequence of the color-singletness requirement since the effective number of internal degrees of freedom 
for a color-singlet QGP is drastically reduced with decreasing volume, as it can clearly be seen from the curves 
and as it has been shown in 0- Thus, the pressure of the QGP phase is lower at a given temperature, and the 
mechanical Gibbs equilibrium between the two phases would then be reached for T C (V) > T c (oo) . 

It can also be noted that, for decreasing volume, while the height of the peak decreases, its width gets 
larger. To see this in more details, we illustrate in the following the second derivative of the order parameter 
(F)(T,V))" = d 2 W^ y)> , which reaches its extrema at temperatures T X {V) and T 2 {V). The width of the 
transition region can simply be defined by the gap between these two temperatures, i. e., 5T(V) — 72(V)— Tj.(V). 
The variations of (f) (T,V))" with temperature for various sizes are illustrated in Fig. (3), from which it can 
clearly be seen that the gap between the two extrema decreases with increasing volume. 

Finally from the obtained results for the finite size response functions, four finite size effects can be observed: 
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• the rounding effect of the discontinuities 

• the smearing effect of the singularities 

• the shifting effect of the transition point 

• the broadening effect of the transition region, 

and to study quantitatively the volume dependence of these effects, a numerical scaling analysis is carried 
out, which will be presented in the next section. 

4. FINITE-SIZE SCALING ANALYSIS 
4.1. Finite-Size Scaling 

In statistical mechanics, it is known that only in the thermodynamic limit are phase transitions characterized 
by the appearance of singularities in some second derivatives of the thermodynamic potential, such as the 
susceptibility and the specific heat. For a first order phase transition, the divergences are (5-function singularities, 
corresponding to the finite discontinuities in the first derivatives of the thermodynamic potential, while for a 
second order phase transition the singularity has a power-law form. For the thermal QCD DPT studied in this 
work, ^-singularities appear in the susceptibility % an d in the specific heat density c at the thermodynamic 
limit. In finite volumes, these (5-functions are found to be smeared out into finite peaks. To the finite size 
effects, four useful characteristic quantities can be associated as illustrated on Fig. (4), which are the maxima 
of the peaks of the susceptibility Xt** (^0 anc ^ the specific heat density c™ ax (V), the shift of the transition 
temperature tt{V) — T c (V) —T c (oo) and the width of the transition region 8T (V) . Each of these quantities can 
be considered as a signature which may anticipate the behavior in the thermodynamic limit, and is expected to 
present a scaling behavior described by a power law of the volume characterized by a Scaling Critical Exponent. 
For a first order phase transition, the set of power laws is: 



x ^ x (V) ~ V~i 
c™ ax (V) - V a 
5T(V) ~ V-° 

t t (V)=T c (V)-T c {oo)~V- 



(21) 



and it has been shown in the FSS theory [Illll2llT3| that in this case, the scaling exponents 9, A, a and 7 are 
all equal to unity, and it is only the dimensionality which controls the finite size effects. 

At a second order phase transition, the correlation length diverges as: £ oc \T — T c \~ u , and then we predict 
the same power law behavior as for a first order transition, but with different scaling critical exponents, usually 
given as: 

X max ry^ _ yl !v 
c max {V) „ ya/v 

ST (V) ~ V~ x l v [ZZ) 
t t (V) = T c (V) - T c (00) ~ V^l v . 

The values of the scaling critical exponents may thengive an indication on the order of a PT and are usually 
used as a criterion for the determination of this latter pj EI • 



4.2. Numerical Determination of the Scaling Critical Exponents for the Thermal DPT 

In the following, we use a FSS analysis to recover the scaling exponents 9, A, a and 7 for the thermally driven 
DPT. For this purpose, we proceed by studying the behavior of the response-function maxima, their rounding as 
well as the shift of the effective transition temperature with varying volume. Let's note that the determination 
of the location of the maxima of the finite size peaks as well as their heights is done in a numerical way, and 
this yields a systematic error, which is estimated and given for the determined scaling exponents. 
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4-2.1. Susceptibility, Specific Heat and Smearing Scaling Exponents 

The data of the maxima of the rounded peaks of the susceptibility | Xt | max 00 an d the specific heat 
density c™ ax (V) are plotted versus volume in Figs. (5-left) and (5-right) respectively, and their linearity with 
V can clearly be noted. A numerical parametrization with the power-law forms: | \t | max 00 ~ ^ 7 an d: 
c max ^y-j ^ ya ^ gj ves ^he values of the susceptibility scaling exponent: 7 = 1.01 ± 0.03, and the specific heat 
scaling exponent: a = 1.007 ± 0.031, where the associated errors are systematic ones. 

Fig. (6) illustrates the plot of the results of the width ST (V) with the inverse of the volume, and their fit to 
the power law form: ST (V) ~ V~ e . The obtained smearing scaling exponent is: 8 = 1.03 ± 0.03. 



4-2.2. The Shift Scaling Exponent 

For the study of the shift of the transition temperature tt(V) — T c (V) — T c (00), we need to locate the 
effective transition temperature in a finite volume T c (V). A way to define T c (V) is as being the location of the 
maxima of the rounded peaks of the susceptibility and the specific heat, shifted away from the true transition 
temperature T c (oo). Results of the shift of the transition temperature obtained in this way are plotted in 
Fig. (7) versus inverse volume. The shift critical exponent obtained from a fit to the form: tt(V) ~ V~ x , is: 
A = 0.876 ±0.041. 



4-2.3. The Shift from the Fourth Binder Cumulant 

Another way for locating T c (V) is to consider the fourth order cumulant of the order parameter proposed in 
[l2| and defined as: 

B 'W^-|||' < 23 » 

which presents a minimum value at an effective transition temperature whose shift from the true transition 
temperature is of order V -1 for a first order transition. This cumulant has been proven to be a suitable 
indicator of the order of the transition in a finite volume, since a nonvanishing value of (| — £? 4 min ) signals a 
first order transition, whereas for a second order transition, lim (| — B4) = even at the transition point. 

Indeed, for a first order transition the quantity lira (3 — -B4) vanishes at all points apart from the transition 

point, and (| - S 4 min ) measures the latent heat [ll H3. ITS IT^ . 

The expression of (T, V) as function of the integral terms, for this case of the deconfmement transition, is 
after calculation: 

B A (T V) = 1 - — - - — - - L ° 5 ~ L ° 4 ^ - 12 ^ 03 ~ ^ L ° 2 Z Lal ^ (24) 

3 (2L13 — 2L03 — 2io2 — ^oi) 

Fig.(8-left) illustrates the variations of the fourth cumulant of the order parameter with temperature for 
various volumes, and shows that the locations of the minima in finite sizes T m ; n (V) are shifted to higher 
values from T c (00). Data of the shift of the transition temperature obtained in this way are plotted in Fig.(8- 
right) versus inverse volume, and the scaling shift critical exponent obtained from a fit to the form: t' t {V) = 
T m in (V) ~ T c (oo) - V~ x ' , is: A' = 0.883 ± 0.043. 



4.2.4. Parametrization of the Order Parameter and the Susceptibility 

Within the model used in this work, the order parameter in the limit of infinite volume being equal to 1 
below the transition temperature and zero above, it can then be expressed in a simple way using the Heaviside 
step-function as: 

(F) (T, V — » 00)) = 1 - 6 (T - T c (00)) . (25) 
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Such a parametrization has been used in [TtI ] , for a similar case of a system of coexisting hadronic matter and 
QGP phases. Using one of the known mathematical representations of the smoothed step function (T — T c ) , 
the order parameter may then be expressed as: 

(i) (T, V)) = Ul- tanh ( T ~ T ;^ } ) ) , (26) 



2 V V r T (v) 

where T c is the effective transition temperature and IV the half-width of the rounded transition region, which 
leads to the susceptibility expression : 

X (T, V) = — ~\ . (27) 

2T T {V) cosh 2 (^p) 

The parametrisation choice (26) is the most accepted physically, and it can be understood phenomenologically 
in the context of the Double Gaussian Peaks model where the obtained expressions of the order parameter and 
the susceptibility are very similar to (26) and (27) 0,0]. 

An illustration of such parametrizations of the order parameter at the volume V = 4000/m 3 , and the 
susceptibility at the volume V = 900/m 3 are presented in Figs.(9-left) and (9-right) respectively, and the 
parameters T c (V) and Ft (V) obtained from each fit are given. 

The results for the width of the transition region ST f lt W (V) and the shift of the transition temperature 
T ^fit(i) ^y-j f rom j-jjg parametrization of the order parameter are plotted with inverse volume in Fig. (10) and 
those from the parametrization of the susceptibility, ST ^**( 2 ) (V) and rj lt ^ (V) in Fig. (11). Their fits to 

power law forms: tJ U{1) (V) ~ V~ Xl , r T flt{2) (V) ~ U~ A2 , ST WO (V) ~ and ST flt ^ (V) ~ V' 02 give: 

Ax = 0.830 ± 0.013 , A 2 = 0.857 ± 0.006 , 9 1 = 0.990 ± 0.015 and 6 2 = 1.032 ± 0.003 . 



5. CONCLUSION 



Our work has shown the influence of the finitcness of the system size on the behavior of some response functions 
in the vicinity of the transition point. The sharp transition observed in the thermodynamical limit, signaled by 
discontinuities in the order parameter and in the energy density at a transition temperature T c (oo), is rounded 
off in finite volumes, and the variations of these thermodynamic quantities are perfectly smooth on the hole 
range of temperature. The delta function singularities appearing in the first derivatives of these discontinuous 
quantities, i.e., in the susceptibility and specific heat density, are then smeared out into finite peaks of widths 
ST (V) . The maxima of these peaks occur at effective transition temperatures T c (V) shifted away from the 
true transition temperature T c (oo). A FSS analysis of the behavior of the maxima of the rounded peaks of the 
susceptibility x™ ax (^0 an d the specific heat density c™ ax (V), the width of the transition region ST (V), and 
the shift of the effective transition temperature relative to the true one tt (V) = T c (V) — T c (oo), shows their 
power-law variations with the volume characterized by the scaling critical exponents 7, a, 9, and A respectively. 
Numerical results for these scaling exponents are obtained and are in good agreement with our analytical 
results: 7 = a = # = A = l obtained in except for the shift critical exponent which slightly deflects 
from the analytical value 1. This may be due to the difficulty of locating accurately the peaks of x> c an d B4, 
especially in the case of large volumes for which the peaks become very sharp. However, a first estimate of the 
scaling critical exponents for the thermal DPT within this model has been obtained, and our result expressing 
the equality between all scaling exponents and the space dimensionality is a consequence of the singularity 
characterizing a first order phase transition and agrees very well with the predictions of other FSS theoretical 
approaches [2(j and with the results of both lattice QCD [2l| and Monte Carlo models calculations. 
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Figure Captions 



Fig. 1: Temperature variation of the order parameter (top) and the energy density normalized by T 4 (bottom) 
at /i = 0, for different system volumes. 

Fig. 2: Plot of (Left) the susceptibility \ (T, V) and (Right) the specific heat density c(T,V), versus tem- 
perature for different system volumes. 

Fig. 3: Plot of the second derivative of the order parameter (f) (T,V)}" = d ^r^^ versus temperature for 
different system volumes. 

Fig. 4: Illustration of the finite size behavior of the susceptibility \ {T, V), the specific heat density c (T, V) 
and the second derivative of the order parameter dx/dT. 

Fig. 5: Variation of (Left) the susceptibility maxima | \t | max (^0 an d (Right) the specific heat density 
maxima cr max (V) with volume. 

Fig. 6: Variation of the width of the temperature region over which the transition is smeared 5T(V) with 
inverse volume. 

Fig. 7: Variation of the shift of the transition temperature tt(V) (from the maxima of \ (T, V) and c (T, V)) 
with inverse volume. 

Fig. 8: (Left) Plot of the fourth Binder cumulant B4 (T, V) versus temperature for different system volumes, 
and (Right) variation of the shift of the transition temperature t' t (V) (from the minimum of B4 (T, V)) with 
inverse volume. 

Fig. 9: Illustration of the fits of the order parameter to the form (26) at the volume V = 4000/to 3 , and of 
the susceptibility to the form (27) at the volume V = 900/m 3 . 

Fig. 10: Variations of (Left) the width of the transition region STf u ^(V), and (Right) the shift of the 
transition temperature t£**^(V), obtained from the parametrization in eq. (26), with inverse volume. 

Fig. 11: Variations of (Left) the width of the transition region ST? 1 *^ (V), and (Right) the shift of the 
transition temperature t£**^(V), obtained from the parametrization in eq. (27), with inverse volume. 
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